Spectroscopy of spin-orbit quantum bits in indium antimonide nanowires 
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Double quantum dot in the few-electron regime is achieved using local gating in an InSb nanowire. 
The spectrum of two-electron eigenstates is investigated using electric dipole spin resonance. Singlet- 
triplet level repulsion caused by spin-orbit interaction is observed. The size and the anisotropy of 
singlet-triplet repulsion are used to determine the magnitude and the orientation of the spin-orbit 
effective field in an InSb nanowire double dot. The obtained results are confirmed using spin blockade 
leakage current anisotropy and transport spectroscopy of individual quantum dots. 
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The spin-orbit interaction (SOI) describes coupling be- 
tween the motion of an electron and its spin. In one di- 
mension, where electrons can move only to the left or to 
the right, the SOI couples this left or right motion to ei- 
ther spin-up or spin-down. An extreme situation occurs 
in what is called a helical liquid [1] where, in the pres- 
ence of magnetic field, all spin-up electrons move to the 
left and all spin-down electrons to the right. As proposed 
recently [2, 3], a helical liquid in proximity to a supercon- 
ductor can generate Major ana fermions [4]. The search 
for Majorana fermions in ID conductors is focused on 
finding the best material in terms of a strong spin-orbit 
interaction and large Lande ^-factors. The latter is re- 
quired for a helical liquid to exist at magnetic fields that 
do not suppress superconductivity. High ^-factors of the 
order 50, strong SOI and the ability to induce supercon- 
ductivity put forward InSb nanowires [5, 6] as a natural 
platform for the realization of ID topological states. 

The SOI can be expressed as an effective magnetic field 
B$o that depends on the electron momentum. An elec- 
tron moving through the wire undergoes spin precession 
around B$o with a 7r rotation over a distance I so called 
the spin-orbit length (see Fig. 1(a)). The length I so is 
a direct measure of the SOI strength: a stronger SOI re- 
sults in a shorter I so- I n this letter, we use spin spectra of 
single electrons in quantum dots [7] to extract I so an d the 
direction of i?so- ^ n quantum dots, the SOI hybridizes 
states with different spin [5, 8, 9]. For a single electron, 
the SOI- hybridized spin- up and spin-down states form 
a spin-orbit qubit [10, 11]. For two electrons SOI hy- 
bridization induces level repulsion between singlet and 
triplet states. The resulting level-repulsion gap between 
the well-defined qubit states can be used to measure the 
SOI: the gap size is determined by l$o [5, 8, 9] and the 
gap anisotropy indicates the direction of B$o [12-14]. 

Double quantum dots in InSb nanowires are defined by 
local gating (Figs. 1(b), 1(c)). A finite voltage is applied 
across the source and drain electrodes; and the current 
through the nanowire is measured. Five gates under- 



neath the wire create the confinement potential and con- 
trol the electron number on the two dots [9, 15]. We focus 
on the (1,1) charge configuration (Fig. 1(d)), in which 
both the left and the right dot contain exactly one elec- 
tron, each of them representing a qubit [10, 11, 16-18]. 

The qubit eigenstates are described by the Kramers 
spin-orbit doublet ff and 4- These two states are su- 
perpositions of spin-up and spin-down, and of several of 
the lowest orbital states [20]. Similar to the case of pure 
spin states, a magnetic field B induces a Zeeman split- 
ting Ez = 9^bB between the Kramers doublets, where 
g is the effective Lande ^-factor for a given direction of 
and jiB is the Bohr magneton. The two qubits in 
the (1,1) configuration can either form a Kramers singlet 
state S(l,l) or one of the three triplets T + (l,l), T (l,l) 
and T_(l,l). The states of the qubits are prepared us- 
ing Pauli spin blockade [10, 11, 17, 18, 21] (Fig. 2(a)), 
which relies on the tunneling process from (1,1) to the 
(0,2) spin singlet S(0,2) (note that T(0,2) state is at 5 
meV above S(0,2) and therefore inaccessible for B = 0). 
When the two electrons form a triplet state, tunneling of 
the left electron to the right dot is prohibited by selection 
rules. This absence of tunneling initializes the qubits in 
the so-called blocked (1,1) state and thereby suppresses 
the current of electrons passing through the double dot. 
Leakage current can occur due to hybridization of T(l,l) 
states with S(0,2) induced by SOI and by spin mixing 
between T(l,l) and S(l,l) due to hyperfine interaction 
[8, 15, 22, 23]. 

Transitions between qubit states are induced by a.c. 
electric fields via electric dipole spin resonance (EDSR) 
[10, 11, 16, 25-27]. Voltages at microwave frequencies are 
applied to the left plunger (LP) gate (Fig. 2(a)). The 
oscillating electric field wiggles the electronic orbits. This 
periodic motion results, via SOI, in a rotation of the spin 
[10, 11]. When the microwave frequency is on resonance 
with the double dot level transitions, EDSR can assist in 
overcoming spin blockade thereby increasing the current 
through the double dot. We map out this current increase 
as a function of microwave frequency / and B (Fig. 2(b)). 
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FIG. 1: (color online) (a), An electron moving with momen- 
tum k through the wire experiences a spin-orbit field .Bso 
which rotates the spin by ty after a distance I so- Vector 
E indicates likely direction of the electric field. In the case 
of spin-orbit coupling due to structural inversion asymmetry, 
Bso oc E x k [19]. (b), Schematic of a double quantum dot in 
an InSb nanowire. Red color indicates regions of the nanowire 
which are not depleted by gates. Gates LB, CB, and RB de- 
fine the left, central and right barriers. Gates LP and RP 
are the left and right plungers used to control the electron 
number on each dot. (c), Scanning electron microscopy of 
a nanowire device similar to the one used in the measure- 
ments, (d), Charge stability diagram of the double dot for 
source-drain voltage V s d= 1 mV. Typical charging energy is 
10 meV. Numbers in brackets correspond to the charge occu- 
pation on the left and the right dots. The inset shows the 
charge stability diagram near the (1,1) — > (0,2) charge transi- 
tion for V s d= 5 mV. The detuning axis e is indicated by the 
dashed arrow. 



For weak interdot tunnel coupling the spectrum is de- 
termined by the energies of individual qubits. At B=0 
all four states are degenerate and non-blocked due to fast 
decay to singlet state induced by hyperfine interaction 
[22]. At finite B, parallel configurations (ft, ft) = T + (l,l) 
and (J|, J|) — T_(l,l) split in energy and become blocked, 
while the other two configurations (J|, ft) and (ft, J|) re- 
main non-blocked. EDSR induces transitions between 
'parallel' and 'anti parallel' configurations, resulting in 
an on-resonance current as observed in Fig. 2(b). The 
slopes of the two "V" shaped resonances determine the 
^-factors of the right and left dots, \qr\ = 29.7±0.2 and 
\gi\ = 32.2±0.2 for this plot. Moreover, the g- factors 
of both dots are highly anisotropic as revealed by the 
EDSR spectroscopy for different field orientations (Fig. 
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FIG. 2: (color online) (a), Left: blocked parallel configuration. 
ip is the angle between nanowire axis nw and B. Right: mi- 
crowaves applied to LP gate (red) induce EDSR. Tunneling 
to the right dot is allowed when the left qubit is rotated to an- 
tiparallel configuration, (b) V-shaped EDSR resonances with 
slopes providing gL and gR for <p =130 ° and V s d = 8 mV. 
Larger ^-factor was assigned to the larger dot, i.e to the dot 
with smaller orbital energy (orbital energy is 5 meV for the 
left dot and 7.5 meV for the right dot). V-shaped lines with 
half the slope are two-photon transitions. Enhanced current 
around B= is due to spin mixing in the absence of mi- 
crowaves (see [24] section S2). Resonances at constant / are 
due to photon-assisted tunneling enhanced by cavity modes. 
(At each frequency the maximum current is normalized to 
lpA and a constant offset is subtracted for clarity.) Inset 
shows energy spectrum of weakly coupled double dots with 
arrows illustrating the observed transitions, (c), Current ver- 
sus / and ip for B = 35 mT. Vertical axis on the left is 
rescaled to g=hf/ ' \ibB. (At each field a constant current off- 
set is subtracted for clarity. ) White arrows over grey cylinders 
indicate B-field orientation with respect to nanowire in panels 
(b) and (c). 
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2(c)). The observed anisotropy is likely determined by 
the details of confinement [28, 29] since the g-factor in 
bulk zincblende InSb is expected to be isotropic. 
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FIG. 3: (color online) (a), Current, in color, versus / and B 
for detuning e ~ 0.5 meV ( V s d =-5 mV). Dashed lines are fits 
to a model described in the [24] section S4. Line colors match 
transitions indicated in panel (c). (At each frequency a cur- 
rent offset is subtracted for clarity.) (b), Diagram illustrating 
a strongly coupled double quantum dot realized by applying 
a more positive voltage to the central gate, (c), Energy dia- 
gram deduced from (a) and used to extract the S-T spin-orbit 
gap Ago . Arrows indicate transitions observed in (a). In the 
absence of coupling, the triplet and the singlet state would 
simply cross as indicated by dashed lines. 

When we increase the interdot tunneling (Fig. 3(b)), 
the (1,1) states hybridize with S(0,2) resulting in level 
repulsion between spectral lines. In the absence of SOI, 
only states with the same spin can hybridize e.g. S(l,l) 
with S(0,2). SOI, however, also enables hybridization 
between the singlets and the triplets [7, 9, 23, 30] (Fig. 3; 
see also Fig. 4(e)). All observed transitions in Fig. 3(a) 
can be identified using a simple model which takes into 
account the hybridization between the (1,1) triplets and 
S(0,2) (see [24] section S4). The four avoided crossings 
observed in Fig. 3(a) correspond to the same double 
dot spin-orbit gap Ago between T_(l,l) and the singlet, 
as illustrated in Fig. 3(c). The quantitative comparison 
with the model allows us to estimate the spin-orbit length 
I so = 230 ± 40 nm (see [24] section S5). 

The observed singlet-triplet gap is highly anisotropic 
(Fig. 4). The gap A^q is largest when B is parallel to 
the nanowire axis nw' ^so snrm ks as the direction of 
B is rotated in the sample plane (Fig. 4(b) and 4(c)). 
Finally for B JL nw the gap disappears (Fig. 4(d)). For 
this orientation the resonance line corresponding to the 
T + (l,l) to singlet transition becomes straight indicating 
the absence of level repulsion between T_(l,l) and singlet. 
In addition, the visibility of the T + (l,l) — >• T_(l,l) tran- 
sition vanishes, suggesting that both T + (l,l) and T_(l,l) 
states are completely blocked for this field orientation. 
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FIG. 4: (color online) a, As the left electron tunnels to the 
right it experiences a field Bso- (b-d) The avoided crossing 
in the EDSR spectrum as in Fig. 3(a) for three directions of 
B: Lp =170°; ip = 110°; and lp = 90° (V sd =-5 mV). (At each 
magnetic field an offset is subtracted for clarity.) (e), Tran- 
sitions between (1,1) states and S(0,2) at finite B. The two 
singlet states are hybridized due to tunnel coupling. T + (l,l) 
and T_(l,l) are coupled to S(0,2) due to Bso- This SOI in- 
duced coupling scales as \Bso x B\ for small B [23]. (f-h), / 
versus e and B for the same orientations of B as in (b-d) with 
microwaves off. (i), Extracted values of Ago (see [24] section 
S6) and / at B = 20 mT and e = 0.5 meV (green dot in panels 
(f-h)) as a function of ip. Solid line is a fit to Ago = Aso|cos 
(cp- (f )\ with A SO =5.2±0.3 fieV and cp = 1°±5° . The error 
bars are determined by the width of EDSR resonance. 



The observed anisotropy of AgQ confirms the spin 
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orbit origin of the singlet-triplet level repulsion (see also 
[24] section S3). The gap AgQ* is expected to be pro- 
portional to \B SO x (B/B)\ [23, 31, 32]. When the two 
fields are aligned, singlet and triplet states cannot mix 
and therefore the spin-orbit gap closes (Fig. 4(d)). From 
the observed anisotropy we conclude that Bso points per- 
pendicular to the nanowire and is parallel to the substrate 
plane (Figs. 4(i) and 4(a)). 

The knowledge of Bso orientation provides a substan- 
tial increase in the fidelity of the initialization and read- 
out of spin-orbit qubits [10]. The fidelity is presently 
limited due to unwanted transitions from T + (l,l) and 
T_(l,l) to S(0,2) induced by SOI. When B and B SO are 
misaligned, T + (l,l) and T_(l,l) are coupled to S(0,2) 
(Fig. 4(e)) [23]. The unwanted transitions are mani- 
fest in the d.c. current through the double dot at finite 
magnetic fields (Figs. 4(f), 4(g), 4(h)) [15, 30]. For an 
ideal readout and initialization no current flows after ei- 
ther T + (l,l) or T_(l,l) state is occupied. When B is 
aligned with I?so? T + (l,l) and T_(l,l) become decou- 
pled from S(0,2) and d.c. current is expected to vanish. 
This dramatic suppression of d.c. current is observed 
for B ±n w (Fig. 4(h)). Importantly, both Agg and I 
show almost identical angle dependence further confirm- 
ing that the singlet-triplet hybridization due to SOI is 
absent when B\\ Bso (Fig. 4(i)). 

Given the direction of Bso we can analyze the ori- 
gin of the spin-orbit interaction in InSb nanowires. The 
field Bso depends on the electron momentum k. In a 
simple physical picture, during the interdot tunneling, 
the momentum k is along the nanowire, which is grown 
in the [111] crystallographic direction. In zincblende 
InSb the spin-orbit interaction has two contributions, the 
bulk-inversion asymmetry term (BIA) and the structure- 
inversion asymmetry term (SIA). However, for k || [111] 
the BIA term is expected to vanish [19], and therefore 
the SIA contribution should dominate. The field Bso 
due to SIA is orthogonal to both the momentum and the 
external electric field (Fig. 1(c)). The electric field is 
likely perpendicular to the substrate since the symmetry 
of confinement in the nanowire is broken by the substrate 
dielectric and voltages on the gates. Therefore the direc- 
tion Bso -L nw and in the substrate plane is consistent 
with the SIA spin-orbit interaction. 

We compare the results obtained from EDSR spec- 
troscopy with the spectrum of (0,2) states (Fig. 5(a)) 
[5, 9, 30]. The SOI hybridization of S(0,2) and T+(0,2) 
states leads to a single dot spin-orbit gap A|q . Since the 
energies of the (0,2) states are too large to be accessed 
with microwaves (A$t ~ 5 meV at 5=0), we use the low- 
est energy T + (l,l) level as a probe of the (0,2) spectrum. 
By changing detuning we move T + (l,l) with respect to 
the (0, 2) levels. When T + (l,l) is aligned with either 
S(0,2) or T + (0,2), an increase in d.c. current is observed 
(Fig. 5(b)) [8]. The level repulsion between T + (0,2) and 



S(0,2) is observed at B w 2T (Fig. 5(c)). The single 
dot gap is also strongly anisotropic reaching the smallest 
value for B JL nw (Figs. 5(d), 5(e) and 5(f)). The spin- 
orbit length I so — 310 ± 50 nm estimated from A|q is 
in agreement with the value obtained using EDSR. 
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FIG. 5: (color online) (a), Two electrons in the right quantum 
dot. The separation of the two electrons in the triplet state is 
of the order of the dot size, (b), Charge stability diagram close 
to (1,1) (0,2) transition at B= 1.4 T, for V sd = 7 mV and 
B\\n w . Transitions T+(l,l) ->> S(0, 2) and T+(l,l) -> T+(0,2) 
are indicated by yellow and gray arrows, (c), Resonances 
corresponding to T+(l,l) -> S(0,2) and T+(l,l) -> T+(0,2) 
as a function of B for (p = lS0°. Colors from dark blue (low) 
to red (high) in panes (c), (d) and (e) indicate values of 61/ ds 
in arbitrary units, (d) and (e), Avoided crossing for (p=180° 
and ip=90°. The dashed lines are fits to the model from Ref. 
[9]. (f), The gap A§§ as a function of (p. Solid line is a fit 
to Aig = Ago cos 2 (ip - (p )cos 2 + sin 2 <9 with A' so = 230 
± 10 /xeV, ip = 2° ± 5° and = 10° ± 3°. The error bars are 
determined by average linewidth corresponding to T+(l,l) 
S(0,2) and T+(l,l) -> T + (0,2) transitions. Note that the 
anisotropy of A§§ depends on the relative positions of the 
two electrons in the right dot which may be different from 
nanowire axis. Out-of-plane Bso angle 6 therefore may be 
non-zero due to confinement details of the right quantum dot. 
Measurements at the (1,1) — >> (2,0) transition yielded the same 
in-plane anisotropy for the left dot (data not shown). 

Recent proposals for experimental detection of Majo- 
rana bound states in hybrid nanowire-superconductor de- 
vices require wires with strong spin-orbit coupling [2, 3]. 
Besides InSb, indium arsenide (InAs) and p-type sili- 
con/germanium (Si/Ge) nanowires [33] are among most 
promising material systems for this purpose. Majo- 
rana states are expected to appear at the boundaries of 
the topological superconducting phase. The topological 
phase is predicted to occur if: (i) E z > A and (ii) E top , 
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A > T. Here A is the superconducting gap, E top is the 
gap of the topological phase and T is the temperature. 
Due to large ^-factors in InSb nanowires first requirement 
is satisfied at low magnetic fields even if large gap super- 
conductors such as niobium are used (A ~ 5K). This is 
a clear advantage since low magnetic fields are preferen- 
tial in order not to suppress superconductivity. The size 
of the topological gap E top ~ 2y/EsoA is determined by 
the bulk SOI splitting E SO = ftV( 2m e*so) W- Here h 
is the Planck constant and m* ~ 0.015 m e is the effective 
electron mass (ra e is the electron mass) . We can estimate 
Eso ~ 0-5 K and E top « 3 K for the case of ballistic one- 
dimensional transport. While Eso is expected to be an 
order of magnitude larger for p-type Si/Ge wires [33] the 
E so ~ 0.1-0.3 K is similar for InAs wurtzite nanowires 
[24] (m* w 0.042-0.06 m e for wurtzite InAs [34]). Note 
however that besides strength of SOI experimental de- 
tails such as quality of semiconductor-superconductor in- 
terface as well as disorder may in the end determine the 
most promising material system. Finally we note that 
the anisotropy measurements (Fig. 4 and 5) suggest the 
orientation B\\nw to be optimal for observing Majorana 
states since the maximum mixing of the SOI-split bands 
occurs for B JL B$o and the superconductivity is sup- 
pressed least when B is in the substrate plane. 

We would like to thank J. Danon, Y. Nazarov, M. Rud- 
ner, D. Loss, F. Hassler and J. van Tilburg for discussions 
and help. We acknowledge help with the measurement 
software from R. Heeres and P. de Groot. This work has 
been supported by ERC, NWO/FOM Netherlands Orga- 
nization for Scientific Research and through the DARPA 
program QUEST. 



[1] Stfeda, P. and Seba, P., Phys. Rev. Lett. 90, 256601 
(2003). 

[2] Lutchyn, R. M., Sau, J. D., and Das Sarma, S., Phys. 

Rev. Lett. 105, 077001 (2010). 
[3] Oreg, Y., Refael, G., and von Oppen, F., Phys. Rev. Lett. 

105, 177002 (2010). 
[4] Kitaev, A. Y., Phys. Usp. 44, 131 (2001). 
[5] Nilsson, H. A., Caroff, P., Thelander, C, Larsson, M., 

Wagner, J. B., Wernersson, L.-E., Samuelson, L., and 

Xu, H. Q., Nano Lett. 9, 3151 (2009). 
[6] Nilsson, H. A., Samuelsson, P., Caroff, P., and Xu, H. Q., 

Nano Lett. 12 228 (2012). 
[7] Schreiber, L. R., Braakman, F., Meunier, T., Calado, V., 

Danon, J., Taylor, J. M., Wegscheider, W., and Vander- 

sypen, L. M. K., Nature Commun. 2, 556 (2011). 
[8] Pfund, A., Shorubalko, I., Ensslin, K., and Leturcq, R., 

Phys. Rev. B 76, 161308 (2007). 
[9] Fasth, C, Fuhrer, A., Samuelson, L., Golovach, V. N., 

and Loss, D., Phys. Rev. Lett. 98, 266801 (2007). 
[10] Nadj-Perge, S., Frolov, S. M., Bakkers, E. P. A. M., and 



Kouwenhoven, L. P., Nature 468, 1084 (2010). 
[11] Nowack, K. C, Koppens, F. H. L., Nazarov, Y. V., and 

Vandersypen, L. M. K., Science 318, 1430 (2007). 
[12] Takahashi, S., Deacon, R. S., Yoshida, K., Oiwa, A., Shi- 

bata, K., Hirakawa, K., Tokura, Y., and Tarucha, S., 

Phys. Rev. Lett. 104, 246801 (2010). 
[13] Kanai, Y., Deacon, R. S., Takahashi, S., Oiwa, A., 

Yoshida, K., Shibata, K., Hirakawa, K., Tokura, Y., and 

Tarucha, S., Nature Nanotech. 6, 511 (2011). 
[14] Golovach, V. N., Khaetskii, A., and Loss, D., Phys. Rev. 

5 77,045328 (2008). 
[15] Nadj-Perge, S., Frolov, S. M., van Tilburg, J. W. W., 

Danon, J., Nazarov, Y. V., Algra, R., Bakkers, E. P. 

A. M., and Kouwenhoven, L. P., Phys. Rev. B 81, 

201305(R) (2010). 
[16] Pioro-Ladriere, M., Obata, T., Tokura, Y., Shin, Y. S., 

Kubo, T., Yoshida, K., Taniyama, T., and Tarucha, S., 

Nature Phys. 4, 776 (2008). 
[17] Petta, J. R., Johnson, A. C, Taylor, J. M., Laird, E. A., 

Yacoby, A., Lukin, M. D., Marcus, C. M., Hanson, M. P., 

and Gossard, A. C., Science 309, 2180 (2005). 
[18] Koppens, F. H. L., Buizert, C., Tielrooij, K. J., Vink, 

I. T., Nowack, K. C., Meunier, T., Kouwenhoven, L. P., 

and Vandersypen, L. M. K., Nature 442, 766 (2006). 
[19] Winkler, R., Spin- Orbit Coupling Effects in Two- 

Dimensional Electron and Hole Systems. Springer, 

Berlin, (2003). 

[20] Flindt, C., Sorensen, A. S., and Flensberg, K., Phys. 

Rev. Lett. 97, 240501 (2006). 
[21] Ono, K., Austing, D. G., Tokura, Y., and Tarucha, S., 

Science 297, 1313 (2002). 
[22] Koppens, F. H. L., Folk, J. A., Elzerman, J. M., Hanson, 

R., van Beveren, L. H. W., Vink, I. T., Tranitz, H. P., 

Wegscheider, W., Kouwenhoven, L. P., and Vandersypen, 

L. M. K., Science 309, 1346 (2005). 
[23] Danon, J. and Nazarov, Y. V., Phys. Rev. B 80, 041301 

(2009). 

[24] Supplementary Material accompanies this paper. It 

can be downloaded at http:/ /link. aps.org/supplemental/ 

10. 1103/PhysRevLett. 108.166801 
[25] Bell, R. L., Phys. Rev. Lett. 9, 52 (1962). 
[26] Golovach, V. N., Borhani, M., and Loss, D., Phys. Rev. 

5 74,165319 (2006). 
[27] Laird, E. A., Barthel, C., Rashba, E. I., Marcus, C. M., 

Hanson, M. P., and Gossard, A. C., Phys. Rev. Lett. 99, 

246601 (2007). 

[28] Pryor, C. E. and Flatte, M. E., Phys. Rev. Lett. 96, 
026804 (2006). 

[29] Schroer, M. D., Petersson, K. D., Jung, M., and Petta, 

J. R., Phys. Rev. Lett. 107, 176811 (2011). 
[30] Pfund, A., Shorubalko, L, Ensslin, K., and Leturcq, R., 

Phys. Rev. Lett. 99, 036801 (2007). 
[31] Nowak, M. P., Szafran, B., Peeters, F. M., Partoens, B., 

and Pasek, W. J., Phys. Rev. B 83, 245324 (2011). 
[32] Stepanenko, D., Rudner, M., Halperin, B. L, and Loss, 

D., Phys. Rev. B 85, 075416 (2012). 
[33] Kloeffel, C., Trif, M., and Loss, D. Phys. Rev. B 84, 

195314 (2011). 

[34] De, A., and Pryor, C. E., Phys. Rev. 5 81, 155210 (2010). 



